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Abstract 



Context. Recent Suzaku X-ray observations of the ejecta-dominated supernova remnant W49B have shown that in the global spectrum 
there is a clear indication for the presence of overionized plasma whose physical origin is still under debate. 

Aims. In order to ascertain the physical origin of such a rapidly cooling plasma, we focus on the study of its spatial localization within 
the X-ray emitting ejecta. 

Methods. We confirm the presence of a saw-edged excess (interpreted as a strong radiative recombination continuum) in the global 
spectrum above 8 keV, emerging above the ionization-equilibrium model. We produce a hardness ratio map to determine where the 
plasma is overionized and we perform a spectral analysis of the regions with and without strong overionization. 
Results. We find that the overionized plasma is localized in the center of the remnant and in its western jet, while it is not detected in 
the bright eastern jet, where the expansion of the ejecta is hampered by their interaction with a dense interstellar cloud. 
Conclusions. The location of overionized plasma suggests that the inner ejecta are rapidly cooling by expansion, unlike the outer 
ejecta, for which expansion is hampered by interstellar clouds seen in H2 
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1. Introduction 

X-ray observations of young supernova remnants (SNRs) allow 
us to probe the physical and chemical conditions of the shock- 
heated ejecta, and to study the physical processes involved in 
their interaction with the circumstellar ambient medium. 

W49B is one of the brightest ejecta-dominated SNRs ob- 
served in X-rays, where it shows a jet-like morphology. Its global 
, X-ray spectrum, characterized by intense emission lines from 

■ He-like and H-like ions of overabundant metals (Si, S, Ar, Ca, 
' Fe, and also Cr and Mn) indicates a plasma at collisiona l ion- 
, izat ion equilibrium (CI E, see, for example, Fuiimoto et aT] | 19951 

■ an d iHwang et a02000h . These results are con firmed by spatially 
' resolved spectroscopy with XMM-Newton (iMiceli et al.l l2006l 

hereafter M06) which reveals a significant Ni overabundance 
(Ni/NiQ - lO^j) in the center. M06 found similar chemical com- 
position and temperatures of the ejecta in the central and eastern 
regions, while in the western part the abundances (in particu- 
' lar for Fe) and the temperatures are lower Both the comparison 
of the observed ejecta abundances with the yields predicted by 
explosive nucleosynthesis models (performed by M06) and the 
multip ole expansion spatial analysis (performed by Lo pez et alj 
I2009ah concur in indicating a core-collapse origin for W49B. 

T he easter n border of the remnant shows bright radio emis- 
sion dLacev e t al. 2001), and is spatially coincident with a 
shocked molecular H2 cloud that hampers the expansion of the 
ejecta in the east direction ( being about three ord er of magnitude 
denser than the ejecta, see iKeohane et al.ir2007l and M06), thus 
distorting the eastern jet southward. Molecular hydrogen and [Fe 



II] and radio emission have been observed in the south-western 
edge of W49B, while in the center the [Fe Ilj emission reveals a 
barrel-Uke st ructure "surrounding" the jet ( Keohane et al.l2007h . 
IMiceli et a n (12008) have shown that W49B seems to be the re- 
sult of an aspherical jet-like supernova explosion with explosion 
energy ^ 1.5 X 10^' erg and mass of the shocked ejecta ~ 6 Mq. 
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iKawasaki et al.l (l2005h analyzed the global spectrum of 
W49B and measured the intensity ratio of the H-like to the He- 
like Ka lines of Ar and Ca, finding that the ionization tempera- 
ture is larger than the electron temperature, Tg, thus claiming 
the presence of overionization in W49B. M06 found T- > 
for Ca and no evidence of overionization for Ar. Nevertheless, 
they also pointed out that the estimate of the electron temper- 
ature derived from the global spectrum is not reliable, given 
that it is significantly lower than that obtained from the spa- 
tially resolved spectral analysis of homogeneous regions. M06 
also performed the same analysis on the uniform central region 
of W49B. Finding no evidence for overionization in the spec- 
trum, they concluded that the temperatures of the ^r and Ca ions 
are consistent with T^. Recentlv. lOzawa et alJ (120091 (hereafter 
O09) analyzed the S uzaku spectrum extracted from the whole 
remnant and revealed the presence of a saw-edged bump above 
8 keV. They demonstrated that the bump is associated with a 
strong radiative recombination continuum of iron and produced 
an accurate spectral model to describe both the recombination 
continuum and the lines and to obtain diagnostics for the ove- 
rionized ejecta (simila r radiative reco mbination features have 
been also observed by Yamaguchi et al. 2009 in IC 443, but for 
Si and S). O09 derived that, for Fe, the ionization temperature is 
kT^ ~ 2.7 keV, while kT,, is ~ 1.5 keV. Moreover, O09 estimated 
that the overionized plasma has the same emission measure as 
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the bremsstrahlung emitting plasma, thus suggesting a common 
origin for the two components. These resuhs strongly indicate 
that the Fe ejecta in W49B are cooling so fast that they are overi- 
onized. Nevertheless, there are still a few important open issues. 

First of all the electron temperature derived from the global 
spectrum by O09 is systematically lower than that observed in 
the uniform spectral regions discussed in M06, where kTe ranges 
between ~ 1.8 keV a nd ~ 3 keV, and in the accurately selected 
regions presented in iLopez et"al] (l2009b) (based on Chandra 
data), where kTi, = 1.8 - 3.7 keV. This discrepancy in kT^ may 
be due to the fact that the global spectrum originates from phys- 
ically non-uniform regions of W49B with different temperatures 
and abundances. Secondly, it is important to verify if all the 
ejecta in W49B are overionized. While M06 did not find ove- 
rionization effects for Ar an d Ca, O09 detected overionization 
for Fe. Notice that lLopez et al. (2009b) have shown that the iron 
morphology is very distinct from that of other elements, being 
more asymmetric and more segregated (and localized in the cen- 
tral and eastern parts of the remnant), while Ar and Ca appear 
well mixed and more isotropic. So, in principle, it is possible 
that the physical conditions in the Fe-rich ejecta are different 
from those of the other ejecta. Finally, a spatial localization of 
the overionized plasma may contribute to ascertain the physical 
origin of its rapid cooling, still not understood. 

Here we present the analysis of archive XMM-Newton obser- 
vations of W49B specifically devoted to constrain the spatial dis- 
tribution of the overionized plasma. The data analysis procedure 
is shown in Sect. |2] and our results are presented and discussed 
in Sect.[3]and Sect.lH respectively. 

2. Data processing 

The XMM-Newton data analyzed here consist of the two obser- 
vations 0084100401 and 0084100501 (PI A. Decourchelle) al- 
ready presented in detail in M06 (see their Sect. 2 and Table 1). 
We re-processed the data using the Science Analysis System and 
following the same procedure and criteria described in M06. The 
images are a superposition of the pn images of the two observa- 
tions and are background-subtracted, vignetting-corrected, and 
adaptively smoothed following the recipe described in M06. In 
the images, the background is subtracted by normalizing blank 
sky observations to the level of flux of the science observation in 
the 10-12 keV (12-14 keV) band for the MOS (PN) cameras. 

For spectral analysis we used the same background event 
files used in M06 and the background spectrum was extracted 
from the same region positions on the CCD and subtracted from 
the source spectrum. The background flux is always a few per- 
cents of the signal in the 8-10 keV band; as an example, for the 
pn camera in the spectral region Y (see SectO, it is 5 4% of the 
total signal in the 8-10 keV band and % 9% of the total flux 
in the 8-12 keV band. We first analyzed the spectra extracted 
from the two observations independently, and verified that they 
are compatible. Then, for each camera, we summed the spec- 
tra extracted from the two observations. Spectra were rebinned 
to achieve a signal-to-noise ratio per bin > Scr and Ancillary 
Response Files were produced with the SAS ARFGEN task. 
The spectral fits were performed simultaneously on the MOSl, 
MOS2, and pn spectra using XSPEC (Arnaud 1996). We adopt 
the best-fit model described in M06. The model consists of two 
absorbed thermal component with oversolar abundances (Si, S, 
Ar, Ca, Fe, and Ni) in ionization equilibrium, plus two narrow 
gaussians to model the Cr and Mn lines. We are interested in the 
hard X-ray emission, so we fit the spectra in the 4.4 - 12 keV 
band. M06 have found small scale variations in the best-fit pa- 



rameters of the cold component and, since here we will focus on 
relatively large regions, we expect that these inhomogeneities 
may be an issue in the spectral fittings. However, in the hard 
4.4 - 12 keV band the cold component contribution is negligi- 
ble. Since it is not possible to constrain the parameters of the 
cold component and the Si, S, Ar, and Ca abundances (whose 
emission lines are all below 4.4 keV), we derive the best-fit tem- 
perature and emission measure of the cold component, the Nh, 
and the Si, S, Ar, and Ca abundances by fitting the spectra in 
the 1-9 keV band and then we freeze them to their best-fit val- 
ues when we perform the fits in the 4.4 - 12 keV band. All the 
reported errors are at 90% confidence. 

3. Results 




channel energy (keV) 

Figure 1. pn (upper) and MOS (lower) global spectra of W49B 
in the 4.4 - 12 keV band, together with their best-fit model and 
residuals. The best-fit model consists of two absorbed thermal 
components with oversolar abundances in collisional ionization 
equilibrium. 

We first extract the spectrum of the entire remnant and we 
model it in the 4.4 - 12 keV band by adopting the best-fit model 
described in Sect. |2] We find a best-fit electron temperature 
kTi = 1.72+° °9 keV, quite similar to that reported by O09. The 
spectrum in Fig. [T] shows quite clearly the presence of a saw- 
edged excess above ~ 8.3 keV, first associated with a radiative 
recombination continuum by O09. We also get a poor fit in the 
Fe (at ~ 6.7 keV) and Fe+Ni (around 8 keV) lines in qualitative 
agreement with O09. Even though the excess is weaker in the 
XMM-Newton spectra than in the S uzaku data (where there are 
several spectral bins at more than 5 sigmas from the model, see 
O09), still it is significant. 

To investigate the spatial distribution of the overionized 
ejecta, we produce a hardness ratio map, i. e. a map of the ratio 
of the pn count-rate in the high energy (8.3 - 12 keV) band to that 
in the bremsstrahlung (4.4-6.2 keV) band (upper panel of Fig.|2] 
H-S map). A high ratio indicates the presence of radiative recom- 
bination continuum and/or a hardening of the bremsstrahlung. 
We remark that if the overionization were not present we would 
expect large values of the H-S map in the eastern end of the 
remnant, where the ejecta temperature is larger (see M06 and 
Lopez et al. 2009b). Instead, in our H-S map the low ratios are 
to the East. To quantify the deviations from CIE accounting for 
the temperature inhomogeneities, we calculated the H-S values 
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Figure 2. Upper panel: vignetting-corrected and background- 
subtracted pn map of the ratio of the count-rate in the hard 
(8.3 - 12 keV) and bremsstrahlung (4.4-6.2 keV) bands (i. e. H- 
S map). The image is adaptively smoothed to a signal-to-noise 
ratio 10 and the bin-size is 12". The white contour lines mark the 
EPIC count-rate (at 33% and 66% of the maximum) in the 1-9 
keV band, and the red contours are at 10%, 3 3%, and 66% of th e 
maximum of the radio image at 327 MHz jLacev etalj|2001h . 
The region selected for the spectral analysis is bounded by the 
black ellipse. The 2 pc scale is obtained assuming a distance of 
8 kpc. Lower panel: ratio of the observed H-S values to that ex- 
pected in the seven spectral regions (whose position is indicated 
in the upper panel) discussed in M06 assuming CIE. 



expected assuming CIE in the seven regions discussed in M06, 
{{H - S)f'^, (i - 1, ...,7), computed using the kT2j reported in 
Table 3 of M06) and compared them to those obtained in our 
H-S map, (H - 5)°**. The lower panel of Fig. |2] shows the ra- 
tios (H - S)fy(H - If the H-S variations were due to 
inhomogeneities in the plasma temperature, these ratios would 
be consistent with one, while values > 1 indicate that the hard 
flux is larger than that expected in equilibrium of ionization. In 
regions I and 2 (at the eastern edge of the remnant) the hard- 
ness ratio is in agreement with that expected in CIE, while in 



regions 3-6 (at the central-western part of W49B) it is signifi- 
cantly larger In region 7 (to the South- West) the difference is 
lower, but still significant. Figure |2] therefore clearly shows that 
the hard radiative recombination continuum is stronger in the 
center and in the western region of W49B, while at the eastern 
edge the H-S values are consistent with CIE. As a cross-check, 
we also produced an "excess image" (not shown here) by sub- 
tracting from the pn image in the 8.3 - 12 keV band the underly- 
ing bremsstrahlung contribution (estimated by fitting the global 
spectrum in the 4.4 - 6.2 keV band) and we obtained the same 
results. 

The image analysis therefore indicates that the ejecta are 
overionized in the center and in the western region of W49B, 
while in the bright eastern edge of the remnant, where the ejecta 
are interacting with a dense molecular cloud, there is no signif- 
icant excess. There are weak indications of overionization also 
in the south-western part of the remnant, characterized by strong 
radio and IR emission and faint X-rays. 

As a further check we extract a spectrum from the black el- 
lipse shown in Fig. |2] where we have localized strong overion- 
ization effects (we call this region "Y") and we compare it with 
the spectrum extracted from the complementary part of the rem- 
nant (i. e. the entire remnant minus the overionized region, here- 
after region "N"). The two spectra are shown in Fig. [3] together 
with their best-fit models and residuals. To estimate the signif- 
icance of the differences in the two spectra we imposed the N 
spectrum to have the same spectral bins as the lower-statistics Y 
spectrum (where the signal-to-noise ratio per bin is > 5cr). The 
best-fit parameters are shown in Table [1] 

Figure[3]clearly confirms that the hard X-ray excess is dom- 
inant above 8.3 keV in region Y, while it is much less significant 
in region N. The last row of Table[T]also confirms that the excess 
is very large only in region Y. We remark that the low signifi- 
cance of the excess in region N is not due to statistics effects, 
since in this region the number of counts in the 4.4 - 12 keV 
band is much larger than in region Y (~ 52000 counts in region 
N vs. ~ 36000 counts in region Y). Both the image and spec- 
tral analysis therefore indicate that the overionization effects are 
mainly localized in the central-western part of the remnant. 



Table 1. Best-fit parameters for the spectra extracted in region 
Y (black ellipse in Fig.|2| and region N (i. e. the entire remnant 
minus region Y). The best-fit model and the fitting procedure are 
described in Sect. |2l All errors are at the 90% confidence level. 



Parameter 


Region Y 


Region N 


kT2 (keV) 
reduced ;^'^/d. o. f. 


1 70+U.UJ 

-^-0.02 

2 29+0.03 
5.8+^1 

^•'-0.7 

1.74/722 


1.53 ±0.02 

43+°'" 

"•^-'-0.01 

4.7 ± 0.3 

9± 1 
1.57/722 


reduced x^/<i. o. f.(8.3 - 12 keV)'*"* 


6.1/38 


2.2/38 



* Emission measure per unit area. 

** Abundances derived in CIE, Discrepancies with O09 are 
due to the CIE model adopted here. 

*** Contribution to the;^'^ of the 8.3 - 12 keV spectral bins. 



4. Discussion 

Our analysis confirms the presence of a radiative recombination 
continuum in the high energy X-ray spectrum of W49B. As ex- 
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Figures. Upper panel: pn (green) and MOS (black and red) 
spectra of region Y (black ellipse in Fig. |2]i together with their 
best-fit model and residuals. Lower panel: same as upper panel 
for region N (the complementary part of the remnant). We im- 
posed the N spectrum to have the same spectral bins as the lower- 
statistics Y spectrum (where the signal-to-noise ratio per bin is 
> 5a-). The best-fit parameters are shown in Table[T] 



plained by O09, this indicates that the ejecta (at least the Fe ions) 
are overionized and are rapidly cooling to reach the ionization 
equilibrium. Nevertheless, we have shown that the bright east- 
ern ejecta of W49B do not present tracers of overionization. It 
is important to stress that the eastern jet of ejecta cannot ex- 
pand freely because a large and dense molecular cloud ha mpers 
the expansion in the east direction (Keohane et al.l 12007 ). The 
western jet, instead, is expanding undisturbed, since it does not 
interact with any dense structure. 

The overionization due to a rapid cooling of the ejecta is 
therefore observed only where the ejecta can expand freely. Such 
rapid cooling is then likely associated with the adiabatic ex- 
pansion of the central-western jet. A similar scenario has been 
invoked by lYamaguchi et al.l (l2009l) to explain the overioniza- 
tion effects in the supernova remnant IC 443. We expect that 
the overionized Fe ejecta experienced a strong heating in the 
early phases of the SNR evolution. This is quite reasonable 
considering that the explosion occurred in a very complex en- 
vironment, likely shaped by winds from the progenitor stars, 
and with the circumstellar material showing a barrel-like mor- 
phology (Keohane et al. 2007). The interaction of the remnant 
with such environment can produce strong reverse and reflected 
shocks that can efficiently heat and ionize the ejecta. The sub- 



sequent rapid expansion of the ejecta in the center and to the 
West then produced their rapid cooling. We have estimated that 
a temperature T ~ 5 x 10^ K and J nedt ~ 2.5 x lO'^ s cm"-' 
would be necessary, in order to have equal ionization time (to 
H-like Fe) and electron heating time. Nevertheless, an accurate 
analysis requires a modeling of the hydrodynamic evolution of 
the plasma including the effects of non-equilibrium of ionization 
(e.g. Reale & Orland o 2008). 

Not all the ejecta of W49B show signatures of overioniza- 
tion for the Fe ions and Fig.|2]shows that in regions with similar 
(and large) Fe abundances, like the eastern jet and the central re- 
gions (see M06 and Lopez et al. 2009b), there are different states 
of ionization (equilibrium to the East and overionization in the 
center), while the central and wes tern regions (that ha ve diff'er- 
ent Fe abundances, see M06 and iLopez et al.ll2009b ) are both 
overionized. We conclude that, while we derive overionization 
effects for the Fe ions, not all the Fe-rich ejecta are overionized. 
On the other hand we remind that in the center of W49B no ove- 
rionization eff'ects are visible for the Ar and Ca ions (M06). We 
expect, therefore, diff'erent ionization temperatures and different 
ratios kT,/kTe for diff'erent species and, at least in the case of 
Fe, for diff'erent regions of the remnant. In principle, it is also 
possible that the temperature of the electrons responsible for the 
Fe recombination is different from that of the electrons respon- 
sible for the bremsstrahlung emission in the 4.4 - 6.2 keV band. 
However, from the global spectrum of W49B, O09 found that 
both the radiative recombination and the bremsstrahlung con- 
tinua can be associated with electrons at the same temperature. 
Nevertheless, a spatially resolved study appears necessary in fu- 
ture works. 

Another important issue is related to the determination of the 
value of the ion and electron temperatures in the ejecta. O09 have 
modeled both the radiative recombination continuum and line 
emission by deriving the ionization temperature from the global 
spectrum, finding kT^ ~ 2.6 keV that is larger than the electron 
temperature kT^, ~ 1.5 keV. In our large regions we also find sim- 
ilar values of kTg, as shown in Table[T] Nevertheless, these values 
are significantly lower than those derived from smaller regions 
by M06. For example, our region Y is approximately the union 
of regions 3, 4, 5, and 6 in M06, where they find kTg ~ 1 .8 - 2.6 
keV that is significantly larger than the value kT,, ~ 1 .75 keV re- 
ported in Table[r| It is then important to analyze spectra extracted 
from relatively small regions in order to constrain accurately 
the ion and electron temperature distribution. Unfortunately, the 
available statistics does not allow us such a detailed analysis and 
we then stress the importance of obtaining deeper observations 
of W49B to address this point. 
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